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20.  ABSTRACT  (Continued) 

only  means  by  which  codes  that  predict  ELF  propagation  in  nuclear  environments 
can  be  compared  with  experiment. 

A full-wave  propagation  code  similar  to  predictive  codes  developed  under 
DNA  auspices  is  used  to  compute  ELF  propagation  anomalies  due  to  Starfish,  PCA 
events,  and  a large  number  of  model  x-ray  flares.  The  results  of  these  calcu- 
lations are  compared— to  the  extent  possible— to  the  available  data.  Despite 
the  inadequacies  of  the  data,  we  are  able  to  verify  key  theoretical  predictions. 

Good  agreement  is  demonstrated  between  code  predictions  and  ELF/VLF  noise 
intensity  during  solar  flares.  Both  experiment  and  theory  consistently  show 
that  the  vast  majority  of  flares  cause  the  intensity  to  increase  at  frequencies 
between  30Hz  and  1 kHz,  decrease  between  5 kHz  and  10  kHz,  and  increase  between 
10  kHz  and  30  kHz.  This  agreement,  which  occurs  over  three  frequency  decades, 
strongly  implies  the  correctness  of  the  predictive  codes. 

Although  the  Schumann-resonance  band  (8-to-30  Hz)  is  somewhat  lower  than 
the  ELF  communications  band  (45-to-80  Hz),  both  bands  propagate  in  the  TEM 
waveguide  mode.  Data  on  Schumann  resonances  may  therefore  be  meaningfully  I 
compared  with  the  outputs  of  ELF  propagation  codes.  The  results  of  such  a 
comparison  made  in  this  report  semiquantitatively  support  the  validity  of  the 
codes.  Specifically,  the  numerical  calculations  correctly  predict  that  either 
Starfish  or  the  PCAs  would  cause  the  peak  frequencies  of  the  Schumann  resonances 
to  shift  downward  by  0.5  Hz  to  1 Hz  and  the  Q's  of  the  resonances  to  decrease. 
The  shift  in  resonant  frequencies  is  due  to  the  depressed  ionosphere  causing 
the  phase  velocity  of  the  signal  to  decrease  and  the  decrease  in  Q is  caused  by 
an  increase  in  the  attenuation  rate  in  the  earth-ionosphere  waveguide. 
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SUMMARY 

An  extensive  literature  survey  has  identified  a number  of  low-frequency 
noise  measurements  relevant  to  the  effects  of  ionospheric  disturbances  on 
ELF  propagation.  These  measurements  pertain  to  the  behavior  of  earth- 
ionosphere  cavity  resonances  and  ELF/VLF  noise  intensity  at  the  time  of 
1)  the  Starfish  high-altitude  nuclear  detonation,  2)  more  than  thirty  solar 
flares,  and  3)  two  polar-cap-absorption  events  (FCAs) . Until  controlled- 
source  ELF  experiments  are  carried  out  under  disturbed  conditions,  these 
collateral  data  provide  the  only  means  by  which  codes  that  predict  ELF 
propagation  in  nuclear  environments  can  be  compared  with  experiment. 

A full-wave  propagation  code  similar  to  predictive  codes  developed 
under  DNA  auspices  is  used  to  compute  ELF  propagation  anomalies  due  to 
Starfish,  FCA  events,  and  a large  number  of  model  x-ray  flares.  The 
results  of  these  calculations  are  compared — to  the  extent  possible — to 
the  available  data.  Despite  the  inadequacies  of  the  data,  we  are  able 
to  verify  key  theoretical  predictions. 

Good  agreement  is  demonstrated  between  code  predictions  and  ELF/VLF 
noise  intensity  during  solar  flares.  Both  experiment  and  theory  consistently 
show  that  the  vast  majority  of  flares  cause  the  Intensity  to  increase  at 
frequencies  between  30  Hz  and  1 kHz,  decrease  between  5 kHz  and  10  kHz, 
and  Increase  between  10  kHz  and  30  kHz.  This  agreement,  which  occurs  over 
three  frequency  decades,  strongly  Implies  the  correctness  of  the  predictive 
codes. 

Although  the  Schumann-resonance  band  (8-to-30  Hz)  is  somewhat  lower 
than  the  ELF  conmiuni cat ions  band  (45-to-80  Hz) , both  bands  propagate  in  the 
TEM  waveguide  mode.  Data  on  Schuaiann  resonances  may  therefore  be  meaningfully 
compared  with  the  outputs  of  ELF  propagation  codes.  The  results  of  such  a 
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comparison  made  In  this  report  semlquantltatlvely  support  the  validity  of 
the  codes.  Specifically,  the  numerical  calculations  correctly  predict  that 
either  Starfish  or  the  PCAs  would  cause  the  peak  frequencies  of  the  Schumann 
resonances  to  shift  downward  by  0.5  Hr  to  1 Hz  and  the  Q's  of  the 
resonances  to  decrease.  The  shift  In  resonant  frequencies  Is  due  to  the, 
depressed  Ionosphere  causing  the  phase  velocity  of  the  signal  to  decrease 
and  the  decrease  In  Q Is  caused  by  an  Increase  In  the  attenuation  rate  In 
the  earth- Ionosphere  waveguide. 


PREFACE 


The  goal  of  this  report  is  to  verify,  insofar  as  is  possible,  codes 
that  predict  ELF  propagation  in  nuclear  environments.  An  extensive  literature 
search  is  made  to  assemble  atmospheric  noise  measurements  relevant  to  ELF 
propagation  during  ionospheric  disturbances.  These  data  are  compared  with 
calculations  made  using  existing  ELF  propagation  codes. 


3 


-j|ii|jij|^U  UWHIIJ 


f 


TABLE  OF  CONTENTS 


SUMMARY  1 

PREFACE  2 

I.  INTRODUCTION  7 

II.  LONG-WAVE  PROPAGATION  CODES  9 

ELF/VLF  PROPAGATION  9 

APPLICABILITY  TO  SCHUMANN -RESONANCE  DATA  12 

III.  SUMMARY  OF  AVAILABLE  DATA 19 

NOISE  INTENSITY  MEASUREMENTS  19 

Effects  of  X-Ray  Flares  19 

Effects  of  the  Starfish  Hlgh-Altltude  Nuclear  Burst  22 

Effects  of  PCA  Events  23 

STRUCTURE  OF  SCHUMANN  RESONANCES  24 

Effects  of  the  Starfish  Hlgh-Altltude  Nuclear  Burst  24 

Effects  of  PCA  Events  29 

Effects  of  X-Ray  Flares  30 

IV.  COMPARISON  BETWEEN  CODE  PREDICTIONS  AND  AVAILABLE  DATA  31 

NOISE  INTENSITY  31 

Effects  of  X-Ray  Flares  31 

Effects  of  the  Starfish  Hlgh-Altltude  Nuclear  Burst  38 

Effects  of  PCA  Events  39 

STRUCTURE  OF  SCHUMANN  RESONANCES  40 

Effects  of  the  Starfish  Hlgh-Altltude  Nuclear  Burst  40 

Effects  of  PCA  Events  43 

Effects  of  X-Ray  Flares  46 

V.  CONCLUSIONS  47 

REFERENCES  49 


4 


LIST  OF  FIGURES 


I . INTRODUCTION 


It  is  frequently  asserted  that  the  ability  of  current  theory  to  predict 
the  performance  of  ELF  communications  systems  in  nuclear  environments  is 

unproven  because  relevant  experimental  data  do  not  exist.  Strictly  speaking,  j 

such  assertions  are  true  because  no  experiment  using  a controlled  ELF  source  ! 

in  a nuclear  environment  has  ever  been  performed,  nor  is  one  likely  to  be 

performed  so  long  as  the  test-ban  treaty  remains  in  force.  Moreover,  no 

data  have  been  obtained  on  transpolar  signals  from  a controlled  ELF  source 

during  a major  polar-cap-absorption  (PCA)  event,  which  would  exhibit  some 

of  the  characteristics  of  a nuclear  environment. 

i 

Accordingly,  this  report  gives  the  results  of  a program  to  use  collateral 
data  from  low-frequency  noise  measurements  to  verify — Insofar  as  it  is 
possible — the  ability  of  current  codes  to  predict  the  effects  of  ionospheric 
disturbances  on  received  ELF  signal  strength.  Obviously,  this  approach  has 
severe  limitations,  and  the  results  cannot  be  expected  to  be  as  satisfactory 
as  would  those  obtained  from  dedicated  experiments  using  controlled  sources. 

Nonetheless,  as  discussed  below,  it  does  prove  possible  to  verify  semlquantl- 
tatively  some  importtnt  predictions  of  existing  codes. 

f 

An  important  aspect  of  the  research  described  in  this  report  is  an 
extensive  literature  survey  to  accumulate  data  relevant  to  the  excitation 
and  propagation  of  ELF  signals  under  disturbed  conditions.  The  most  useful  | 

I 

of  the  data  reviewed  are  measurements  of  ELF/VLF  atmospheric  noise  intensity  1 

and  the  structure  of  the  earth-cavity  (Schumann)  resonances  during  the  Starfish 
hlgh-altltude  nuclear  test,  several  PCA  events,  and  numerous  sudden  ionospheric 
disturbances  (SIDs).  We  do  not  consider  in  detail  data  from  the  Wisconsin  Test 
Facility  ELF  transmitter,  which  have  hfen  adequately  addressed  by  other  inves- 
tigators (c.g.,  Imhoff,  et  at.,  1976). 
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Section  II  reviews  the  salient  aspects  of  ELF  propagation  theory  and 
codes  based  thereon  and  establishes  the  applicability  of  these  codes  to 
Schtimann  resonance  phenomena;  Sec.  Ill  summarizes  the  data  accumulated 
during  the  literature  search;  Sec.  IV  pr'^sents  computational  results  and 
compares  them  with  the  data  given  In  Sec.  Ill;  Sec.  V summarizes  the 


conclusions. 


II.  LONG-WAVE  PROPAGATION  CODES 


This  section  briefly  summarizes  some  important  aspects  of  currently 
used  long-wave  propagation  codes  and  establishes  the  applicability  of 
these  codes  to  the  Interpretation  of  measurements  of  ELF/VLF  atmospheric 
noise  and  the  Schumann  resonances. 

ELF/VLF  PROPAGATION 

The  detailed  equations  governing  ELF/VLF  propagation  in  the  earth- 
ionosphere  waveguide  are  given  in  a number  of  references  (e.g.,  BuddeK, 
Field,  1970;  Pappert  and  Moler,  1974;  Field,  Lewinetein,  and  Done, 
1970'),  and  will  not  be  repeated  here.  These  references,  among  others, 
also  describe  the  numerical  procedures  employed  by  current  ELF/VLF  codes 
to  solve  these  equations.  These  codes  account  for  the  vertical  inhomo- 
genelty  of  the  ionosphere,  the  horizontal  inhomogeneity  of  the  ionosphere 
along  the  great-circle  connecting  transmitter  and  receiver,  and  the  curva 
ture 'of  the  earth;  but  neglect  horizontal  inhomogeneities  transverse  to 
the  great-circle  path.  Field  {1978)  considers  the  effects  of  transverse 
Inhomogeneltles  and  gives  a correction  term  to  be  used  if  the  transverse 
dimension  of  an  ionospheric  disturbance  is  smaller  than  a Fresnel  zone. 

To  define  symbols  and  illustrate  key  dependences,  the  equation  for  the 
spatial  dependence  of  ELF/VLF  fields  is  given  and  discussed  below. 

Usually,  somewhat  dif f erei\^3quations  are  used  for  ELF  and  VLF  propa 
gat  ion  because  the  WTF  ELF  antenna  is  horizontally  oriented,  whereas  most 
operational  VLF  transmitters  are  vertically  oriented.  Here,  however,  we 
are  concerned  with  atmospheric  noise  generated  by  the  vertical  component 
of  lightning  strokes,  and  a single  equation  may  therefore  be  used.  In 
general,  only  one  waveguide  mode — the  TEM  mode — can  propagate  at  ELF; 
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but  several  modes  can  contribute  to  VLF  fields.  In  this  report,  we  restrict 
attention  to  long  paths,  so  that  higher  modes  are  heavily  attenuated,  and 
only  a single  VLF  mode — the  TM— 1 mode — need  be  considered. 

Subject  to  the  above  conditions,  the  following  equation  describes  the 
vertical  electric  field  at  both  ELF  and  VLF: 


- 120irle 


-wl/A 


✓Xd 


d/a 


sin  d/a 


53/2^^-Bd/(8.7xlO°)  ^-(2irl/X)  (c/v) 


(1) 


In  Eq.  (1),  IL  Is  the  effective  electric  dipole  moment  of  the  source,  X Is 
the  free-space  wavelength,  d Is  the  distance  from  source  to  receiver,  a Is 
the  earth's  radius,  and  c Is  the  vacuum  speed  of  light.  MKS  units  are  used 
except  where  expressly  Indicated  otherwise. 

The  role  of  propagation  codes  Is  to  determine  the  dependence  of  the 
signal,  E,  on  the  state  of  the  Ionosphere.  This  dependence  occurs  solely 
through  two  quantities:  1)  the  excitation  factor.  A,  which  governs  the 
efficiency  with  which  the  waveguide  mode  Is  generated;  2)  the  quantity,  S, 
which  Is  essentially  the  eigenvalue  of  the  relevant  waveguide  mode.  The 
losses  suffered  by  a signal  as  It  propagates  In  the  earth- Ionosphere  wave- 
guide are  given  by  the  attenuation  rate,  B,  which  Is  expressed  In  decibels 
of  loss  per  megameter  of  propagation  (dB/Mm) . The  phase  velocity  of  the 
signal  Is  denoted  by  v.  Both  B end  v are  related  to  S by  the  following 
simple  relations: 

S = + 1 , (2) 

B - 0.18  f Sj,  dB/»ta  , (3) 

c/v  - Sj^  , (4) 

where  f Is  the  wave  frequency. 
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The  quantity  actually  computed  by  long-wave  propagation  codes  Is  the 
so-called  elgencoslne,  C,  of  the  relevant  waveguide  mode.  The  closely  related 
quantity,  S,  Is  then  obtained  from 


This  distinction  between  C and  S Is  Important  In  the  verification  of  codes 
because  and  are  the  quantities  that  can  be  determined  experimentally, 
whereas  C Is  the  quantity  actually  computed.  The  main  point — especially  at 
ELF,  where  neither  the  real  nor  the  Imaginary  part  of  C Is  particularly 
small — Is  that  depends  strongly  on  both  Cj^  and  Thus,  the  correctness 

of  a calculated  value  of  S„  (phase  velocity)  Implies  a correct  value  of  S. 

R 1 

(attenuation  rate),  and  vice-versa. 

The  above  discussion  shows  that  output  of  a long-wave  propagation  code 
Is  essentially  two  quantltl'^:,  A and  S,  and  that  verification  should  Involve 
comparison  of  measured  value  of  these  quantities  with  calculated  results. 
This  procedure  Is  vastly  more  complicated  than  Indicated  by  the  apparently 
simple  equations.  The  calculation  of  S requires  a highly  sophisticated 
Iterative  numerical  technique  that  uses  detailed  profiles  of  Ionospheric 
charged-particle  densities  and  collision  frequencies  as  Inputs.  Moreover, 
nature  Is  seldom  so  cooperative  as  to  provide  a uniform  propagation  path 
characterized  by  a single  value  of  S at  all  locations.  Thus,  even  If 
transverse  Ionospheric  variations  can  be  Ignored,  measured  ELF  data  pertain 
to  a propagation  path  along  which  various  Ionospheric  profiles  and  ground 
conductivities  are  encountered. 

Finally,  we  note  that,  although  all  equations  given  In  this  section 
are  formally  valid  at  both  ELF  and  VLF,  Implicit  differences  exist  between 
the  two  bands.  For  frequencies  below  about  1 kHz,  the  eigenvalue  S and 
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excitation  factor  A correspond  to  the  TEM  mode,  which  has  a somewhat 
different  structure  than  the  TM-1  mode,  to  which  S and  A correspond 
above  a few  kilohertz.  However,  the  numerical  procedures  used  to  calculate 
S and  A are  identical  for  these  two  modes. 

APPLICABILITY  TO  SCHUMANN-RESONANCE  DATA 

The  equations  given  above  are  in  a form  applicable  to  ELF  communica- 
tions frequencies  of  many  tens  of  hertz,  for  which  the  wavelength  is  much 
smaller  than  the  circumference  of  the  earth.  Some  of  the  most  useful  of 
the  available  data  were  measured  in  the  so-called  Schumann  resonance  band. 
This  band  Includes  frequencies  between  about  8 Hz  and  30  Hz,  which  are 
somewhat  lower  than  frequencies  usually  associated  with  ELF  communications. 
Electromagentlc  signals  at  8-to-30  Hz  propagate  in  the  same  TEM  mode  as  do 
45-to-80-Hz  ELF  communications  signals,  and  essentially  the  same  theoretical 
treatment  applies  equally  well  to  the  two  bands.  Thus,  a validation  of 
code  predictions  in  the  Schumann  band  strongly  Implies  that  predictions  in 
the  usual  ELF  band  are  also  correct.  However,  because  wavelengths  in  the 
Schumann  band  can  be  comparable  to  the  earth's  circumference,  it  is  con- 
venient to  write  the  equation  for  the  electric  field  in  a form  that  differs 
somewhat  from  Eq.  (1). 

Galejs  (1972)  shows  that  the  vertical  electric  field  is  given  by 


„ -v  IIL  (kaS)^ 
8,fc  ha2  cos(nkaS) 


P^(-  cos  d/a) 


v/m 


(6) 


where  P is  the  Legendre  function,  k is  the  free-space  wave  number,  h is  the 
nominal  reflection  height  of  the  ionosphere,  and  the  approximation 
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V 


(7a) 


1 


*,Vo 


0.25  + (kaS)  - ^ 


^ kaS  - >s  (7b) 

has  been  made.  Equation  (6)  Is  actually  more  accurate  at  all  frequencies 
than  Eq.  (1),  which  Is  an  approximate  form  convenient  for  use  at  frequencies 
higher  than  about  30  Hz.  Specifically,  by  using  the  large-order  asymptotic 
expansion  of  P^,  neglecting  waves  that  travel  around  the  earth  more  than 
once,  noting  that  A ^ 1/h,  and  performing  some  algebraic  manipulation.  It 
Is  straightforward  (but  tedious)  to  show  that  Eq.  (6)  may  be  approximated 
by  Eq.  (1) — provided  that  the  wavelength  Is  much  smaller  than  the  earth's 
c ircumf  erence . 

Equation  (6)  shows  that  the  spectrum  of  background  noise  would  be 
expected  to  exhibit  peaks  near  frequencies  where  the  denominator  exhibits 
minima;  viz.,  spectral  peaks  should  occur  where 


cos(iTkaS)  = minimum 


(8) 


Equation  (8)  will  be  satisfied  at  a series  of  n complex  frequencies  given 
by  (n  » 1,  2,  ) 


fn  = f + If,  - 

n in  2^5 


Hz 


(9) 


Thus,  the  center  frequency  of  the  n^''  resonant  peak  Is  given  by 


f Hz 

n 2¥  a 


and  the  "Q"  of  the  n*”^  peak  Is  given  by 


(10) 


13 


I 


I 


I 


1^ 


I 


I 

! 


i 


Q 


n 


(11) 


2 2 

In  Eqs.  (10)  and  (11),  use  has  been  made  of  the  fact  that  S.  Recall 

1 K 

that  Is  a measure  of  the  width  of  a resonant  peak  and  becomes  smaller 

If  the  peak  becomes  wider;  l.e.,  if  the  peak  Is  "smeared."  If  the  earth 

and  Ionosphere  were  sharply  bounded  perfect  conductors,  S would  equal 

K 

unity,  and  Eq.  (10)  would  predict  resonances  at  10.6  Hz,  18.3  Hz,  25.9  Hz, 
etc.  Calculations  using  realistic  ionospheres,  however,  show  that  is 
somewhat  larger  than  unity,  reducing  the  predicted  resonant  frequencies  to 
around  8 Hz,  14  Hz,  20  Hz,  26  Hz,  etc.,  under  ambient  conditions.  Figure  1 
shows  a nominal  measured  noise-power  spectrum,  similar  to  that  first 
measured  by  Balser  and  Wagner  (1960),  The  theoretically  predicted  peaks 
are  evident. 

Schumann  resonances  are  particularly  convenient  for  interpretation 
because  f^  and  depend  solely  on  and  S^,  which  essentially  are  the 
phase  velocity  and  attenuation  rate  of  the  signal  (see  Eqs.  (3)  and  (4)). 

Thus,  unlike  measurements  of  noise  intensity  that  depend  on  source  strength 
and  excitation  factor,  measured  values  of  f^  depend  to  a good  approximation 
solely  on  the  phase  velocity,  and  thus  present  no  ambiguity  in  interpretation. 

The  goal  of  this  report  is  to  establish — Insofar  as  is  possible — that 
currently  used  codes  correctly  predict  ELF  signal  changes  associated  with 
ionospheric  disturbances.  Specifically,  a change  in  ionospheric  structure 
will  cause  A and  S to  assume  new  post-disturbance  values,  denoted  by  A + AA 
and  S + AS.  In  a strong  nuclear  environment,  AS^  is  so  large  as  to  create 
a severe  Increase  in  the  attenuation  rate  8,  thus  causing  blackout.  For 
the  disturbances  of  natural  origin  considered  in  Secs.  Ill  and  IV  below. 
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Fig.  1— Typical  measured  power  spectrum  of  atmospheric  noise  in  the 
O-to-40-Hz  band.  Note  the  resonances  at  around  8 Hz,  14  Hz, 
20  Hz,  etc. 


It  Is  valid  to  assume  AS<<S;  and  it  follows  from  Eqs.  (10)  and  (11)  that 


= 1 


I 


f 


Af  = 


n 


and 


(12) 


. (13) 

where  Eq.  (3)  has  been  used  to  relate  S to  B,  and  S-  and  B are  evaluated 

X K 

at  f . Equations  (12)  and  (13)  show  that:  a)  an  Increase  in  S (equivalent 
n R 

to  a decrease  in  phase  velocity)  will  cause  a downward  shift  In  the  resonant 
frequencies;  b)  an  Increase  in  attenuation  rate,  B,  causes  a reduction  In  Q 
(equivalent  to  broadening  of  the  spectral  peaks) . 

The  resonance  condition,  Eq.  (8),  and  the  resulting  Eqs.  (12)  and  (13) 
are  applicable  only  for  a stratified  Ionosphere  where  the  disturbance  is 
such  that  ASj^  and  AB  do  not  depend  on  geographic  location.  Thus,  these 
equations  must  be  generalized  to  be  applicable  to  realistic  disturbances 
of  limited  geographic  extent.  To  accomplish  this  generalization,  we  expand 
the  denominator  of  Eq.  (6)  In  the  well-known  series 


cosdkaS)  " exp  KirkaS  - ,r/2)  . (14) 

m=0 

Note  that  wa  is  the  distance  a wave  travels  In  propagating  halfway  around 
the  earth,  and  irkaS  is  the  phase  change  corresponding  to  this  distance.  It 
is  therefore  evident  from  Eq.  (14)  that  Schumann  resonances  occur  from  the 
constructive  self -interference  of  waves  that  travel  around  the  earth  a 
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large  number  of  times,  while  undergoing  a ir/2  phase  change  each  time  the 
antipode  is  passed.  The  resonance  condition,  (Eq.  (8)),  simply  defines  the 
discrete  frequencies  for  which  the  phase  term  irkaS  assures  that  a wave  will 
constructively  self-lnterfere  after  having  made  an  Integer  number  of  round 
trips. 

The  best  available  data  on  the  effects  of  ionospheric  disturbances  on 
Schumann  resonances  were  taken  during  the  Starfish  hlgh-altltude  nuclear 
test  and  several  PCA  events.  The  resulting  disturbances — although  not 
worldwide — were  far  from  being  localized.  Under  these  conditions.  It  Is 
reasonable  to  use  the  phase-memory  approximation,  and  note  that  the  phase, 
4,  corresponding  to  propagation  from  source  to  antipode  can  be  written 


♦ « irka 


S + 


II 

^ / 


d0AS(e) 


J 


radians 


(15) 


where  6 Is  the  usual  angular  spherical  coordinate  and  the  Integral  Is  taken 
along  the  great-circle  path  containing  source  and  receiver.  Thus,  the 
effective  value  of  the  perturbation  on  S Is 


<AS 


T 

>-/ 


d0AS(e) 


(16) 


where  <AS>  Is  the  average  value  of  AS  taken  over  the  path.  By  substituting 
S + <AS>  for  S In  the  resonance  condition,  Eq.  (8),  It  Is  easily  shown  that 
Eqs.  (12)  and  (13)  become 


Af 


<ASr> 


Sr  " 


(17) 
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and 


<AB> 

6 


(18) 


If  good  data  were  available  on  the  location  of  sources  and  lateral 
Ionospheric  structure  at  the  time  of  the  measured  disturbances.  It  would 
be  reasonable  to  perform  full-fledged  calculations  of  the  average  values 
<AS->  and  <AB>  Indicated  In  Eqs.  (17)  and  (18).  Unfortunately,  the  data 

K 

summarized  In  Sec.  Ill  are  much  too  coarse  to  warrant  a complex  calculation, 
and  there  Is  no  choice  but  to  use  rather  crude  estimates  for  <AS->  and  <AB>. 

K 

Specifically,  we  use 


<ASj^>  FAS^ 


<AB>  ^ FAB 


where  AS^^  and  AB  are  taken  to  be  characteristic  values  along  the  path 
containing  source  and  receiver  and  extending  the  circumference  of  the 
earth,  and  F denotes  the  estimated  fraction  of  this  path  over  which  these 
characteristic  values  obtain.  The  final,  very  approximate,  form  of  the 
equations  used  to  relate  Ionospheric  disturbances  to  changes  In  Schumann 
resonances  thus  becomes 


Af 

n 


FAS„ 


(19) 


AQ 


n 


(20) 
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III.  SUMMARY  OF  AVAILABLE  DATA 


This  section  summarizes  the  results  of  a literature  survey  to  assemble 
data  pertinent  to  the  effects  of  ionospheric  disturbances  on  ELF  signal 
strength.  As  Indicated  in  Sec.  II,  the  behavior  of  Schumann  resonances 
during  ionospheric  disturbances  is  relevant  to  the  verification  of  ELF 
predictive  codes,  although  Schumann  resonance  frequencies  are  lower  than 
those  under  consideration  for  ELF  communications.  Therefore,  this  section 
gives  Schumann  resonance  data  as  well  as  data  for  the  behavior  of  ELF  noise 

intensity  during  various  types  of  disturbances.  In  some  situations,  VLF  i 

noise  data  were  taken  in  conjunction  with  ELF  noise  data,  thus  affording 
an  opportunity  to  interpret  the  relative  behavior  of  both  bands.  Thus, 
although  this  report  is  concerned  with  verification  of  theoretical  predic- 
tions at  ELF,  some  VLF  data  are  given,  where  appropriate. 

NOISE  INTENSITY  MEASUREMENTS 

1 

Effects  of  X-Ray  Flares 

Scientists  from  Japan  have  carried  out  two  measurement  programs  to 
determine  the  change  in  intensity  of  background  ELF/VLF  noise  caused  by 
solar  x-ray  flares.  The  first  set  of  measurements,  made  in  May  1963, 
monitored  ELF  frequencies  of  30  Hz,  100  Hz,  260  Hz,  and  600  Hz;  and  VLF 
frequencies  of  10  kHz,  21  kHz,  and  27  kHz  (Sao  and  Jindoh,  1966).  The 
results  of  those  measurements  are  shown  in  Table  1 , which  presents  the 
change  in  noise  intensity  that  occurred  at  the  onset  of  four  flares.  Note 
that,  in  almost  every  Instance,  the  noise  intensity  increased  somewhat, 
motivating  the  term  SEA  (sudden  enhancement  of  atmospherics)  to  describe 
the  phenomenon. 
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The  second  set  of  measurements,  made  during  1968  and  1969,  was  much  more 
extensive  than  the  first  {Sao,  et  al^ , 1970).  More  than  20  events  were  moni- 
tored, and  more  frequencies  (100  Hz,  170  Hz,  260  Hz,  400  Hz,  570  Hz,  800  Hz 
at  ELF;  5 kHz,  10  kHz,  21  kHz,  27  kHz  at  VLF)  were  included.  Figure  2 summa- 
rizes the  results  of  this  second  set  of  measurements. 

The  experimental  results  shown  In  Table  1 and  Fig.  2 cover  nearly 
30  events  over  a period  of  several  years  and  reveal  a remarkably  repeatable 
phenomenon.  In  every  instance,  the  flares  caused  enhancement  of  noise  in  the 
ELF  band  (800  Hz  and  below)  and  the  upper  VLF  band  (21  kHz  and  27  kHz). 
Conversely,  in  the  lower  VT,F  band  the  noise  was  decreased,  with  the  decrease 
being  much  more  severe  at  5 kHz  than  10  kHz.  This  decrease  in  noise  gave 
rise  to  the  term  SDA  (sudden  decrease  of  atmospherics),  with  the  overall 
effect  of  flares  on  background  noise  denoted  by  SEA/SDA.  In  several 
Instances,  the  SEA/SDA  effects  were  substantial,  exceeding  10  dB  in 
magnitude.  The  SEA/SDA  effect  at  VLF  had  been  noted  by  several  authors 
{Pierce,  1961;  Chilton,  et  al. , 1964;  Crombie,  1966;  Burgees  and  Jones,  1967), 
but  the  SEA  effects  at  ELF — which  is  of  primary  interest  here — appears  to 
have  been  measured  solely  by  Sao  and  his  collegues. 

Table  1 

MEASURED  SEA/SDA  RATIO  IN  DB 
(A  POSITIVE  SIGN  DENOTES  AN  SEA;  A NEGATIVE,  AN  SDA) 

{Sao  and  Jindoh,  1966) 


Date 

Frequency 

30  Hz 

100  Hz 

260  Hz 

600  Hz 

10  kHz 

21  kHz 

27  kHz 

2 May  1963 

+1.7 

+2.1 

+5.9 

+3.8 

No  obser- 
vation 

+1.0 

+3.5 

16  May  1963 

+1  • 6 

+2.0 

+2.7 

+2.9 

-0.9 

0 

+1.5 

20  May  1963 

+1.0 

+1.6 

+1.9 

+3.2 

0 

0 

+1.0 

21  May  1963 

0 

0 

+2.5 

-1.1 

0 

0 

+3.0 

frequency 


Fig.  2--Measured  frequency  spectrum  of  SDA/SEA  {Sao,  et  al.,  1970). 
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Unfortunately,  no  data  are  available  on  the  intensity  or  spectrums 
of  the  Incident  x-rays  responsible  for  the  event  shown  in  Table  1 or 
Fig.  2.  Thus,  it  is  not  possible  to  make  a detailed  calculation  relating 
a specific  flare  to  a specific  SEA/SDA  event.  This  lack  of  detailed  informa- 
tion is  not  really  too  serious,  however,  because  the  remarkable  consistency  of 
the  SEA/SDA  data  shows  that  the  general  features  of  the  phenomenon  are 
virtually  independent  of  the  structure  of  the  flare.  Specifically, 

SEA  was  consistently  observed  below  800  Hz  and  above  20  kHz,  and 
SDA  was  consistently  observed  at  5 kHz  and  10  kHz — despite  the  fact 
that  the  group  of  nearly  30  flares  must  have  included  some  that 
differed  substantially  from  others  in  intensity  and  spectral  character- 
istics. 

Effects  of  the  Starfish  Hlgh-Altltude  Nuclear  Burst 

As  indicated  above,  unambiguous  interpretation  of  the  effect  of  a 
nuclear  burst  on  noise  intensity  is  difficult  because  the  weapon  can 
affect  source  strength  and  the  propagation  characteristics  of  the  earth- 
ionosphere  cavity.  Even  if  data  were  taken  over  a long  enough  time  period 
to  avoid  obvious  transients,  such  as  the  EMP,  one  cannot  be  certain  whether 
to  attribute  a change  in  noise  intensity  to  a bomb-induced  change  in  wave- 
guide structure  or  to  some  long-term  bomb-induced  change  in  source  charac- 
teristics; e.g.,  trapped-partlcle  populations  or  lightning-flash  distribu- 
tions. Nonetheless,  for  completeness,  it  appears  worthwhile  to  summarize 
ELF-noise  intensity  data  taken  at  the  time  of  Starfish.  Starfish  had  a 
yield  of  1.4  HT  and  was  detonated  400  km  over  Johnston  Island  in  the  Pacific 
at  about  0900  UT  on  9 July  1962. 
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Teply,  et  al.  (,1963)  measured  the  ELF  background-noise  intensity  at 
frequencies  between  5 Hz  and  1.5  kHz  at  Tongatapu,  which  is  about  3000-mi 
south  of  ground-zero.  The  data  given  in  Table  2 were  measured  during  time 
periods  subsequent  to  the  first  few  seconds  after  the  burst  and  therefore 
avoid  the  transient  EMP.  In  some  instances,  the  nature  of  the  data  are 
such  (e.g.,  sonograms  or  traces  without  scales)  that  only  qualitative 
Information  can  be  given  in  Table  2. 


Table  2 

SUMMARY  OF  ELF  DATA  MEASURED  AT  TONGATAPU 
BY  TEPLY,  ET  AL.  AT  TIME  OF  STARFISH 


Effects  of  PCA  Events 

Larsen  (1974)  presents  data — taken  at  Troms^,  Norway — on  the  behavior 
of  ELF  atmospheric  noise  during  two  PCA  events.  The  first  set  of  measure- 
ments— taken  continuously  during  the  10-day  period  from  18  October  1969  to 
6 November  1969— monitored  frequency  bands  centered  at  8 Hz,  35  Hz,  75  Hz, 
and  115  Hz.  The  measured  field  strength  exhibited  a A-to-8  dB  diurnal 
variation  that  qualitatively  follows  statistics  for  worldwide  thunderstorm 
activity.  A major  PCA  event  started  at  around  1100  UT  on  2 November  1969 
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and  persisted  for  at  least  five  days.  Despite  the  strength  and  duration 
of  this  event,  no  significant  change  was  noted  In  the  Intensity  and  tem- 
poral behavior  of  the  noise  at  the  four  fre<;uencles  monitored. 

Surprisingly,  a small  PCA  that  occurred  on  24  January  1969  apparently’ 
caused  a significant  change  In  ELF  noise  Intensity,  whereas  the  major 
2 November  event  discussed  above  did  not.  Figure  3 shows  the  measured 
power  spectral  density  between  a few  hertz  and  400  Hz  before,  during,  and 
after  the  24  January  PCA.  The  PCA  appears  to  have  enhanced  the  noise 
somewhat  at  frequencies  below  about  100  Hz,  and  reduced  the  noise  by 
several  dB  at  frequencies  between  100  Hz  and  400  Hz.  The  difference  In 
noise  during  the  PCA  (24  January)  from  that  proceeding  the  PCA  (20  January) 
cannot  be  attributed  to  diurnal  variation  of  thunderstorm  activity,  because 
both  measurements  were  taken  at  about  the  same  time  of  day.  However, 

Larsen  states — and  we  agree — that  the  single  set  of  measurements  shown 
In  Fig.  3 Is  Inadequate  to  show  conclusively  that  the  24  January  PCA  was 
the  cause  of  the  altered  noise  levels. 

STRUCTURE  OF  SCHUMANN  RESONANCES 

Effects  of  the  Starfish  Hlgh-Altltude  Nuclear  Burst 

The  spectrums  of  electromagnetic  background  noise  In  the  5-to-30-Hz 
band  were  measured  at  the  time  of  the  Starfish  hlgh-altltude  nuclear  burst 
by  Balser  and  Wagner  {196 S) , Gendrln  and  Stefant  {1962,  1964),  and  Teply 
et  al.  {196S) , By  comparing  the  locations,  f^,  of  the  spectral  peaks 
before  and  after  the  burst,  the  shift,  Af^,  In  each  peak  caused  by  the 
burst  can  be  estimated. 

The  main  results  of  these  measurements  are  summarized  in  Table  3. 

The  time  constants  of  the  measuring  equipment  varied  from  tens  of  seconds 
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-Power  spectral  density  function  recorded  at  Lavangsdulen,  Norway  in  1969  {Laraen,  1974). 


lY  OF  MEASURED  SHIFTS  IN  SCHUMANN-RESONANCE  FREQUENCIES  CAUSED  BY  THE  STARFISH  BURST 
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(continued) 


Table  3 (continued) 


Authors  state  that  determination  of  spectral  peak  was  "crude. 


to  several  minutes;  and,  to  within  these  timing  uncertainties,  the  Indicated 
changes  In  the  resonant  frequencies  appeared  to  occur  Instantly.  Teply,  et  al. 


did  not  specify  uncertainties  In  their  measured  values  of  but  stated 
that  their  determination  of  f^  subsequent  to  the  shot  was  "very  crude" 
because  of  smearing  at  the  resonant  peaks.  Note  that  the  measurements  were 
made  at  three  very  widely  spaced  locations;  viz.,  Tongatapu,  about  3000-ffll 
south  of  the  burst;  Massachusetts,  about  6000-mi  northeast  of  the  burst; 
France,  more  than  8000  ml  from  the  burst  via  a transpolar  great-circle  path. 

The  data  consistently  show  that  the  Starfish  burst  caused  the  first 


few  resonant  frequencies  to  diminish  by  significant  amounts.  The  magnitudes 
of  the  decreases  depend  on  the  resonance  number,  being  about  1/2  Hz  for  the 
lowest  resonance  and  around  1 Hz  for  the  fourth  resonance.  The  2.4-Hz 
decrease  in  the  third  resonant  frequency,  f^,  reported  by  Teply,  et  al.  is 
much  larger  than  any  of  the  other  reported  frequency  shifts.  In  view  of 
Teply 's  own  statements  regarding  the  imprecision  of  his  spectral-peak  deter- 
mination, we  treat  his  measurement  of  Af^  as  spurious. 

Unfortunately,  no  specific  measurements  were  reported  on  the  effect  of 
Starfish  on  the  Q's  of  the  Schumann  resonances.  Visual  inspection  of  spec- 

trums  presented  by  Teply,  et  al.,  however,  shows  that  the  spectral  peaks  were 

* 

slightly  smeared  by  the  burst.  Also,  the  fifth  and  sixth  resonances. 


which  were  visible  on  the  pre-shot  spectrums,  were  essentially  undefined 
on  the  post-shot  spectrums.  Moreover,  Balser  and  Wagner  (1963)  comment 
that  the  burst  noticeably  degraded  the  quality  of  the  first  few  resonant 
peaks.  Thus,  although  no  quantitative  conclusions  can  be  drawn  regarding 
Q,  we  feel  justified  in  drawing  the  qualitative  conclusion  that  Starfish 
caused  the  Qs  to  decrease  somewhat . 


Our  visual  interpretation  of  these  spectrums  is  not  consistent  with  that 
of  Teply,  et  al.,  who  state  that  the  burst  produced  "no  noticeable  change  in 
the  Q of  the  earth-ionosphere  cavity." 
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Gendrln  and  Stefant  (1964)  present  maps  showing  the  locations  of 
thunderstorm  activity  during  the  months  of  July.  Statistically,  the 
sources  contributing  to  natural  background  electromagnetic  noise  at  the 
time  of  Starfish  would  be  expected  to  be  heavily  concentrated  in  Southeast 
Asia.  Of  course,  the  possibility  exists  that  the  Starfish  burst  altered 
the  geographic  distribution  of  noise  sources. 

Effects  of  PGA  Events 

Nelson  (1967)  presents  data  for  the  frequency  and  Q of  the  lowest 
Schumann  resonance  during  three  PCA  events.  The  measured  changes  in  fj^ 
and  Qj^  are  given  in  Table  4,  along  with  the  dates  and  times  of  the  three 
events.  This  table  also  shows  the  nominal  centers  of  the  thunderstorm 
activity  causing  the  atmospheric  noise  at  the  time  of  the  PCA  events.  As 
would  be  expected,  the  PCA-induced  changes  are  very  similar  to  those 
described  above  for  the  Starfish  burst.  Specifically,  the  resonant  fre- 
quency decreased  by  several  tenths  of  a hertz  and  Q decreased  by  an 
ill-defined — but  noticeable — amount . 

Table  4 

PCA-INDUCED  CHANGES  IN  FREQUENCY  AND  Q OF  LOWEST  SCHUMANN 
RESONANCE;  MEASURED  IN  MASSACHUSETTS 
(Nelson,  1967) 


f^(Pre-PCA) 

f^(PCA) 

Afj 

AQ 

Date 

Local 

Standard 

Time 

Thunderstorm 

Center 

Not  given 

Not  given 

-0.2  Hz 

"slight 

decrease" 

6 July  1966 

2030 

Scattered 

8.4  Hz 

8.0  Hz 

-0.4  Hz 

"moderate 

Increase" 

28  Jan  1967 

1000 

Africa 

8.2  Hz 

7.4  Hz 

-0.8  Hz 

"slight 
decrease  ' 

21  Feb  1967 

1400 

South 

America 

J 


Effects  Of  X-Ray  Flares 


A large  number  of  measurements  have  been  made  by  scientists  from  Japan 
during  SIDs  (Sao^  et  al.,  1971,  1973;  Ogaua,  et  al.,  1966;  OgcaXL  and  Tamaka, 
1970).  Specifically,  the  correlations  of  the  frequencies  and  Qs  of  the 
resonances  with  sunspot  number,  lA-to-SA  x-ray  intensity,  and  magnetic 
Ap-index  were  determined.  For  the  most  part,  these  correlations  are 
inconclusive,  although  an  imaginative  reader  might  discern  an  Inverse 
relationship  between  1)  resonant  frequency  and  sunspot  number  and  2)  Q and 
Ap-lndex.  Nelson  (.1967)  gives  data  for  the  frequency  and  Q of  the  lowest 
resonance  during  a number  of  SIDs  and  concludes  that  individual  events  fail 
to  produce  discernible  effects. 
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IV.  COMPARISON  BETWEEN  CODE  PREDICTIONS  AND  AVAILABLE  DATA 

This  section  compares  the  results  of  theoretical  calculations  with 
the  experimental  data  described  In  Sec.  III.  The  calculations  are  performed 
using  Paclf Ic-Slerra  Research's  full-wave  propagation  code  (see,  e.g.. 

Field,  1969;  Field,  1970;  Field,  Lewinstein,  and  Dove,  1976),  which  Is 
virtually  Identical  with  ELF  predictive  codes  developed  and  exercised 
under  DNA  auspices  (see,  e.g.,  Pappert  and  Moler,  1970),  For  ease  of 
comparison,  the  various  categories  of  data  are  discussed  below  In  the  same 
order  as  In  Sec.  III.  This  ordering,  however,  does  not  Indicate  the  rela- 
tive Importance  of  the  data  categories.  The  data  on  the  effects  of  x-ray 
flares  on  ELF  noise  Intensity,  and  on  the  effects  of  the  Starfish  burst 
and  PCA  events  on  Schumann-resonance  structure,  are  the  most  conclusive 
and  strongly  confirm  key  predictions  of  the  propagation  codes.  On  the 
other  hand,  data  on  the  effects  of  Starfish  and  PCAs  on  ELF  noise  intensity, 
and  on  the  effects  of  flares  on  Schumann  resonances  tend  to  be  Inconclusive, 
and  neither  strongly  support  nor  strongly  disagree  with  code  predictions. 

NOISE  INTENSITY 
Effects  of  X-Ray  Flares 

To  compare  code  predictions  with  the  experimental  results  given  In 
Sec.  Ill,  calculations  of  ELF/VLF  attenuation  rate,  phase  velocity,  and 
excitation  factor  are  made  for  model  Ionospheres  based  on  ambient  conditions 
and  a number  of  assumed  solar-flare  conditions.  The  ambient  daytime  model 
ionosphere  Is  taken  from  the  DNA  Blackout  Handbook  (Knapp  and  Sohuartz,  197SJ; 
the  disturbed  ionospheres  are  calculated  from  assumed  x-ray  flares  and  spec- 
trums  using  lumped-parameter  de-lonolzatlon  equations  and  reaction  rates 
also  given  In  the  Blackout  Handbook.  An  infinite  ground  conductivity, 
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corresponding  approximately  to  a sea  water  path.  Is  assumed.  Emphasis  Is 
placed  on  calculating  flare-induced  ahangee  In  propagation  parameters  rather 
than  on  absolute  values,  because  the  data  given  In  Table  1 (p.  20)  and  Fig.  2 
(p.  21)  denote  flare- Induced  changes  In  field  strength. 

The  SEA  and  SDA  Indicate  that  the  frequency  de,'endence  of  flare-induced 
changes  In  ELF/VLF  noise  Intensity  Is  not  sensitive  to  the  detailed  structure 
of  a given  flare.  Moreover,  these  data  are  not  presented  In  a form  that  per- 
mits direct  comparison  between  a specific  flare  and  a specific  set  of  noise 
measurements.  It  Is  therefore  not  worthwhile  to  use  detailed  models  of 
specific  flares  to  make  calculations.  The  most  profitable  approach  Is  to 
perform  calculations  for  a number  of  nominal — but  realistic — flare  models, 
and  compare  the  ensemble  of  results  thus  obtained  with  the  ensemble  of  data 
given  in,  e.g..  Fig.  2. 

It  Is  not  altogether  unfortunate  that  the  available  data  must  be 
subjected  to  a pseudostatistlcal — rather  than  one-for-one — interpretation, 
because  well-defined  models  of  flares  do  not  exist.  The  x-ray  flux  and 
spectrum  Is  different  for  different  flares,  and  changes  during  a given 
flare.  Specifically,  the  flux  tends  to  Increase  and  the  spectrum  tends  to 

O 

harden  as  a flare  progresses.  X-rays  in  the  l-to-8-A  spectral  band  are  most 
effective  In  Influencing  long-wave  propagation,  because  they  cause  Ionization 
in  the  D and  lower  E regions.  Papers  by  De  Jager  (.1984,  1965),  Friedman  (1963), 
Pounds  (1966),  and  Van  Allen  (1969),  among  others,  summarize  data  on  x-ray 

e o 

flux  and  spectrum  In  this  band.  Integrated  fluxes  between  lA  and  8A  varying 

2 

from  hundredths  to  tenths  of  an  erg/cm  -sec  were  measured.  A wide  variety 

e 

of  spectra  were  noted.  Including  sharp  peaks  in  the  S-to-lO-A  range,  and 
quasi-thermal  spectra  at  blackbody  temperatures  up  to  around  lO^K. 


To  represent  the  above  conditions  in  the  l-to-8-A  band,  we  have  performed 

calculations  for  several  flare  models  having  assumed  Integrated  fluxes  ranging 
2 2 

from  0.05  erg/cm  -sec  to  0.2  erg/cm  -sec,  and  characterized  by  both  monoener- 
getlc  and  blackbody  spectra.  The  monoenergetlc  spectra  were  characterized  by 

O O 

x-ray  energy  peaks  between  1.5  KeV  (8A)  and  3 KeV  (4A) , and  the  blackbody 
spectra  by  "temperatures"  between  5xl0^”K  and  10^ ®K. 

Figure  4 shows  the  attenuation  rate  In  decibels  per  megameter  of 
propagation  calculated  for  the  ambient  daytime  model  ionosphere,  and 
Ionospheres  corresponding  to  model  flares  described  above.  Only  the  envelope 
of  calculated  values  for  the  flare  models  is  shown,  because,  as  discussed 
above,  comparison  of  the  data  with  specific  flares  is  not  possible.  The 
results  given  In  Fig.  4 show  that  the  attenuation  rates  calculated  for  the 
model  flares  are  about  equal  to  ambient  below  100  Hz,  lower  than  ambient 
between  100  Hz  to  1 kHz,  and  between  10  kHz  to  30  kHz,  and  higher  than 
ambient  between  5 kHz  and  10  kHz.  Since  Increased  attenuation  corresponds 
to  signal  degradation,  the  results  in  Fig.  4 are  consistent  with  the  SEA/SDA 
data  in  Table  1 and  Fig.  2. 

To  better  illustrate  the  agreement  between  calculations  and  experiment, 
the  results  from  Fig.  4 have  been  plotted  in  a different  format  In  Fig.  5, 
which  shows  the  differenae  between  attenuation  rates  calculated  for  the  model 
flares  and  those  calculated  for  the  model  ambient  Ionosphere.  The  sign  of 
this  difference  has  been  reversed  to  account  for  the  fact  that  an  Increase 
In  attenuation  rate  corresponds  to  a decrease  in  electromagentlc  field  inten- 
sity. Thus,  Fig.  5 shows  the  flare- induced  enhancement  or  degradation  pet 
megameter  of  propagation  path  exposed  to  the  flare. 

In  addition  to  enhancements  or  degradations  that  are  due  to  propagation 
and  are  proportional  to  exposed  pathlength,  solar  flares  can  affect  the 
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model  solar- flare  environments. 


:hanqes  in  ELF/VLF  attenuation  rate  relative  to  ambient  for  several  model 
environments 


excitation  factor  of  the  dominant  waveguide  modes.  Figure  6 shows  the 
calculated  excitation  factors  for  the  model  solar  flares.  These  exclta- 


i. 


1 
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tlon  factors  are  expressed  In  decibles  relative  to  the  ambient  value.  As 
Is  typical  of  almost  all  types  of  Ionospheric  disturbances,  the  excitation 
factors  are  Increased  over  their  ambient  value.  Note  that  the  enhance- 
ments associated  with  excitation  factors  are  independent  of  pathlength,  as 
opposed  to  the  propagation  anomalies  shown  in  Fig.  5,  which  must  be  multi- 
plied by  total  pathlength  to  obtain  the  change  in  measured  field  strength. 

The  agreement  between  the  computational  results  of  Figs.  S and  6 with 
the  data  of  Table  1 and  Fig.  2 is  striking.  Below  100  Hz,  the  attenuation 
rate  (Fig.  5)  is  essentially  unchanged,  but  the  excitation  (Fig.  6)  is 
enhanced  by  0 to  1 dB.  These  results  are  in  excellent  agreement  with  the 
data  of  Table  1,  which  shows  O-to-2  dB  enhancements  at  30  Hz  and  100  Hz. 

At  frequencies  between  100  Hz  and  30  kHz,  the  structure  of  the  calculated 
results  (Fig.  5)  is  essentially  identical  to  that  of  the  data  shown  in 
Fig.  2;  viz.,  enhancement  between  100  Hz  and  1 kHz,  degradation  between 
5 kHz  and  10  kHz,  enhancement  above  10  kHz.  Moreover,  the  spectral  shape 
of  the  calculated  results  is  very  similar  to  that  of  the  data,  with  the 
strongest  degradation  occurring  at  frequencies  of  a few  kilohertz;  i.e., 
the  strongest  degradation  occurs  just  above  the  cutoff  of  the  lowest  TM 
mode.  By  multiplying  the  calculated  propagation  anomalies  (Fig.  5)  by 
reasonable  disturbed  pathlengths — say,  a few  megameters — and  adding  the 
excitation  anomalies  (Fig.  5),  it  is  easy  to  show  that  the  calculations 
are  in  reasonable  quantitative  agreement  with  the  data  from  Table  2. 

The  data  of  Table  2 contain  no  phase  information,  because  of  the  non- 
controlled  nature  of  the  source.  However,  many  observations  of  signals 
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For  the  above  reason,  the  best  that  can  be  done  is  to  give  a tenta- 
tive qualitative  explanation  of  the  data.  Table  2 shows  that  subsequent 
to  Starfish  the  noise  intensity  increased  at  frequencies  below  about 
30  Hz  and  was  essentially  unchanged  at  higher  frequencies.  Many  existing 
ELF  calculations  (e.g..  Field,  1969)  show  that  the  effect  of  a large  dis- 
turbance such  as  a nuclear  burst  or  strong  PCA  is  to  increase  both  the 
excitation  factor  and  attenuation  rate.  The  Increase  in  excitation  is 
relatively  constant — typically  several  dB — over  the  entire  band,  whereas 
the  increase  in  attenuation  rate  becomes  more  severe  as  the  ELF  frequency 
is  increased.  The  result  is  the  well-known  competition  at  ELF  between 
improved  excitation  and  degraded  propagation  (increased  attenuation)  that 
can  give  either  a net  increase  or  decrease  in  measured  field  strength. 
Qualitatively,  the  improved  excitation  would  be  expected  to  dominate  the 
degraded  propagation  at  the  lower  ELF  frequencies,  and  to  become  less 
Important  relative  to  the  degraded  propagation  at  the  higher  frequencies. 
This  general  behavior  is  consistent  with  the  data  given  in  Table  2. 

Effects  of  PCA  Events 

As  indicated  in  Sec.  Ill,  the  available  data  on  the  effects  of  PCA 
events  on  ELF  noise  intensity  are  Inadequate  for  quantitative  interpreta- 
tion. The  net  effect  of  an  ionospheric  disturbance  on  ELF  noise  Intensity 
depends  strongly  on  the  location  and  orientation  of  the  propagation  path 
relative  to  the  disturbed  region.  Neither  the  geographic  extent  of  the 
PCA  nor  the  location  of  the  noise  sources  relevant  to  the  data  in  Fig.  3 
are  known  with  any  precision.  Thus,  it  is  not  possible  to  calculate  the 
net  effect  of  the  above-mentioned  competition  between  Improved  excitation 
and  degraded  propagation.  Although  surprising  at  first  glance,  the  fact 
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that  a small  PCA  apparently  caused  a greater  change  in  noise  strength  than  did 
a large  event  Is  not  necessarily  at  variance  with  code  predictions.  The  large 
PCA  simply  could  have  caused  an  Increase  in  excitation  that  compensated  for  a 
degradation  in  propagation.  Thus,  we  conclude  that  these  data  are  inconclusive. 

STRUCTURE  OF  SCHUMANN  RESONANCES 

Effects  of  the  Starfish  High-Altitude  Burst 

We  use  Eqs.  (19)  and  (20)  to  compare  theoretical  predictions  of  the  peak 
frequencies,  f^,  and  Q's  of  the  Schumann  resonances  with  data  given  in  Sec.  III. 
Empahsls  is  placed  on  bomb- Induced  shifts,  1^  peak  frequencies  shown 

in  Table  3 (pp.  20-21)  because  the  main  goal  of  propagation  codes  is  to  predict 
signal  changes  caused  by  changes  in  ionospheric  structure.  Interpretation 
of  Schumann-resonance  data  is  quite  unambiguous  because  the  structure  of  the 
resonances  depends  mainly  on  the  state  of  the  ionosphere,  being  Independent 
of  source  strength. 

To  use  Eq.  (19),  we  need  theoretical  calculations  of  the  characteristic 
change,  AS^,  in  relative  phase  velocity  caused  by  the  Starfish  burst.  In 
addition,  we  need  an  estimate  of  the  fraction,  F,  of  the  earth's  circumference — 
drawn  through  source  and  receiver — that  is  subjected  to  levels  of  disturbance 
characterized  by  As  discussed  in  Sec.  II,  ASj^  exhibited  strong  geographic 

disturbance  after  Starfish,  and  a well-defined  value  does  not  exist.  Moreover, 
although  the  main  sources  of  atmospheric  noise  at  the  time  of  Starfish  were 
probably  centered  somewhere  in  Southeast  Asia,  an  accurate  source  location 
is  not  known.  These  facts,  coupled  with  the  diffuse  boundaries  of  the  bomb- 
disturbed  region,  make  an  accurate  value  of  the  fraction,  F,  impossible. 

Thus,  the  best  that  can  be  hoped  for  is  semlquantitative  agreement  between 
code  predictions  and  available  data. 

The  first  step  in  estimating  AS.  is  to  determine  the  value  of  that 

R R 

obtained  prior  to  the  burst.  Even  this  step  is  imprecise  because  the 
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from  VLK  transmitters  in  the  16-to-20  kHz  range  have  been  made  during 
flares  (e.g.,  Cpombie,  1965;  Chilton,  et  at.,  1964;  Buvgese  and  Jones, 

1967;  Pieroe,  1961).  '•’Iiese  measurements  show  that  the  flares  generally 
caused  the  amplitude  to  Increase — in  agreement  with  the  calculations  and 
noise  data  from  Fig.  2 — and  the  phase  to  advance.  Although  not  explicitly 
plotted  here,  phase  calculations  were  made  in  conjunction  with  the  attenua- 
tion rate  calculations  shown  in  Figs.  4 and  5.  In  all  Instances,  the 
calculated  phase  velocity  increased  between  16  kHz  and  20  kHz,  which 
agrees  with  the  phase  advances  shown  by  the  data. 

In  summary,  the  code  predictions  agree  well  with  noise-amplitude 
measurements  made  during  solar  flares  over  the  three  frequency  decades 
from  30  Hz  to  30  kHz.  Even  the  crossover  frequencies  separating  bands  of 
enhancement  from  bands  of  degradation  are  accurately  predicted.  Moreover, 
VLF  phase  calculations  agree  with  measurements  of  flare-induced  phase 
anomalies  between  16  kHz  and  20  kHz.  This  agreement  strongly  supports  the 
validity  of  ELF/VLF  propagation  codes  designed  to  predict  signal  strength 
(and  phase)  changes  caused  by  ionospheric  disturbances. 

j Effects  of  the  Starfish  High-Altitude  Nuclear  Burst 

! The  ELF  noise-amplitude  data  taken  subsequent  to  Starfish  cannot  be 

unambiguously  Interpreted  using  propagation  codes.  This  interpretive  dif- 
ficulty occurs  because  the  codes  deal  only  with  changes  in  attenuation 
rate,  phase  velocity,  and  excitation  factor  caused  by  alterations  in  the 
I structure  of  the  earth- ionosphere  waveguide.  In  addition  to  these  factors, 

^ the  noise  intensity  depends  on  the  location  and  strength  of  noise  sources, 

which  might  well  have  been  altered  by  the  burst. 
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propagation  paths  pertinent  to  Schumann  resonances  completely  encircle  the 
earth,  and  some  averaging  between  day  and  night  conditions  is  needed.  More- 
over, on  the  nighttime  portion  of  the  path,  depends  on  propagation  direc- 
tion. The  average  ambient  values  used  for  taken  from  Galejs  {1972)  are 
as  follows:  1.35  at  8 Hz  (fj),  1.3  at  lA  Hz  1.27  at  20  Hz  (f^), 

and  1.25  at  26.5  Hz  (f^). 

The  full-wave  propagation  code  is  used  to  calculate  values  of  Sj^ 
subsequent  to  the  Starfish  burst.  These  calculated  values  depend  on 
distance  from  the  burst,  and  a single  value  cannot  be  given.  Fortunately, 
once  a relatively  strong  level  of  ionospheric  disturbance  has  been  reached, 
the  Impact  on  ELF  propagation  tends  to  be  saturate;  and  further  Increases 
in  the  level  of  disturbance  cause  relatively  small  changes  in  ELF  propaga- 
tion parameters.  Thus,  it  is  possible  to  define  values  of  S that  corres- 
pond  roughly  to  a fairly  wide  range  of  strongly  disturbed  conditions.  By 
performing  a number  of  numerical  calculations  for  nuclear  environments 
of  various  intensities,  we  find  the  following  values  to  reasonably 
characterize  in  strongly  disturbed  regions:  1.7  at  8 Hz  (fj^);  1.6  at 
14  Hz  (fj);  1.5  at  20  Hz  (f j) ; 1.45  at  26.5  Hz. 

By  Inserting  the  above  disturbed  and  undisturbed  values  for  S in 
Eq.  (19),  it  is  easily  shown  that  the  theoretically  predicted  Starfish- 
induced  shifts  in  resonant  frequencies  become 


Af  'w2F  Hz 
n 

-wSF  Hz 
'W3.5F  Hz 
'w4F  Hz 


for  fj(8  Hz) 
for  fjCK  Hz) 
for  13(20  Hz) 
for  f^(26.5  Hz) 


(21) 


Again,  recall  that  F is  the  fraction  of  the  appropriate  round-the-world 
propagation  path  significantly  affected  by  the  weapon.  Note  that  Starfish 
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affected  enoraous  regions  due  to  line-of-sight  radiation  from  the  burst  and 

subsequent  far-flung  debris,  conjugate  zone  affects,  and  neutron-decay  B 

j 

• particles  Impinging  on  the  Ionosphere  well  beyond  llne-of-slght . Thus, 

. although  It  Is  Impossible  to  accurately  determine  F,  values  In  the  range 

^ of  0.15  to  0.3  seem  reasonable  for  so  widespread  a disturbance  as  Starfish. 

The  larger  values  In  this  range  apply  to  Teply's  Tongatapu  Station,  which 
was  squarely  In  the  conjugate  zone.  Insertion  of  these  values  for  F Into 
Eq.  (21)  gives  good  semlquantltatlve  agreement  with  the  data  In  Table  2 


(p.  23).  Specifically,  the  codes  correctly  predict  that  1)  the  Starfish 
burst  should  have  caused  downward  shifts  In  the  Schumann  resonance  frequen- 


cies, 2)  that  the  magnitude  of  these  downward  shifts  Is  between  several 

* 


tenths  of  a hertz  and  around  a hertz,  and  3)  the  greatest  frequency  shifts 
occur  for  the  higher  order  resonances. 


The  same  numerical  calculations  used  to  obtain  the  values  of  S given 
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above  show  that  a strong  disturbance  such  as  Starfish  will  cause  the  attenua- 
tion rate,  6,  to  Increase  by  a few  tens  of  percent  at  8 Hz  and  by  up  to  a 
factor  of  two  at  26.5  Hz.  Insertion  of  these  values  for  A6/B  into  Eq.  (20) 
show  that  the  codes  predict  a decrease  In  the  Q's  of  the  resonances,  with 
the  greatest  predicted  decreases  occur ing  for  the  higher -order  resonances. 
This  result  Is  consistent  with  the  qualitative  data  Interpretation  given  In 
Sec.  Ill;  viz..  Starfish  apparently  caused  the  Q's  to  decrease  "somewhat." 
Absence  of  better  data  on  the  behavior  of  the  Q's  makes  more  quantitative 
conclusions  regarding  agreement  between  experiment  and  theory  Impossible. 

As  pointed  out  In  the  derivation  given  In  Sec.  II,  however,  the  fact  that 
the  codes  correctly  predict  ASj^  virtually  guarantees  the  correctness  of 
predicted  changes  in  attenuation  rate. 


Recall  that  Teply's  measurement  of  -2.4  Hz  for  Af^  is  probably  spurious. 
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Effects  of  PCA  Events 

Table  2 (p.  23)  summarizes  the  data  on  the  effects  of  PCA  events  on 
the  lowest  Schumann  resonance.  Of  the  three  events  shown,  only  those  of 
28  January  and  21  February  1967  are  considered  here.  The  6 July  event  Is 
difficult  to  Interpret  because  only  partial  data  are  available  and,  because 
of  the  time  of  the  event,  no  reasonable  assumptions  can  be  made  regarding  the 
location  of  major  thunderstorm  centers.  No  Information  was  given  on  the 
structure  of  the  PCAs  or  the  Ionosphere  at  the  time  of  the  measurements 
summarized  In  Table  2.  Thus,  as  was  the  situation  for  the  ELF  noise 
measurements  taken  during  flares,  there  is  no  point  in  attempting  detailed 
correlations  between  specific  flares  and  the  measured  values  of  and  AQ. 
Again,  the  best  that  can  be  done  Is  to  compare  the  data  with  calculations 
based  on  several  nominal  PCA  models. 

Table  5 summarizes  the  results  of  calculations  based  on  three  nominal 
PCA-dlsturbed  daytime  ionospheres  designated  genetically  as  strong,  moderate, 
and  weak.  All  models  are  based  on  actual  events,  albeit  not  the  ones  of 
28  January  and  21  February  1967,  and  are  similar  to  those  discussed  by  Field 
(1970),  The  calculations  are  performed  for  a frequency  of  8 Hz,  which  cor- 
responds closely  to  the  center  frequency  of  the  lowest  Schumann  resonance. 

As  has  been  done  throughout  this  report,  we  concentrate  on  the  changes, 

AS^  and  AS,  produced  by  the  disturbance  rather  than  on  the  absolute  values 
of  Sjj  and  6.  In  Table  5,  ASj^  and  AS  denote  the  difference  between  the  dis- 
turbed and  the  ambient  values.  Table  5 also  Includes  nominal  ambient  values 
for  S„  and  S,  averaged  over  daytime  and  nighttime  conditions. 

K 

As  was  the  situation  for  Interpreting  the  Starfish  data,  it  Is  necessary 
to  estimate  the  factor  F before  applying  Eqs.  (19)  and  (20).  For  the 


43 


21  February  event,  the  main  thunderstorm  centers  were  probably  in  South 
America,  and  the  great-clrcl'e  path  connecting  this  source  location  to  the 
receiver  location  in  Massachusetts  is  essentially  in  the  North-South  direc- 
tion. Thus,  the  round-the-world  path  containing  source  and  receiver  passes 
squarely  through  both  polar  paps,  and — depending,  of  course,  on  the  details 
of  the  particular  PCA — near^  one-third  of  this  path  can  be  assumed  covered 
by  the  disturbance.  Thus,  we  use  F a 0.3  for  the  21  February  PCA  event. 

Inserting  F - 0.3  and  ASj^/Sj^  from  Table  5 into  Eq.  (19)  gives 

♦ 

Afj^  ~ -0.6  Hz  for  moderate-tp-strong  PCAs  and  Af^^  ~ -0.35  Hz  for  weak 
PCAs.  Given  the  many  uncertainties  In  the  calculation,  the  agreement  be- 
tween these  results  and  the  experimental  value  of  Af^  • -0.8  Hz  is  satis- 
factory. Inserting  F « 0.3 ‘and  A6/6  from  Table  5 Into  Eq.  (20)  indicates 
that  the  21  February  PCA  would  have  been  expected  to  decrease  by  be- 
tween 15  percent  and  20  percent.  These  predictions  are  In  qualitative 
agreement  with  the  observed  "slight  decrease"  Indicated  in  Table  2. 

At  the  time  of  the  28  January  1967  PCA,  the  most  likely  location  of 
the  important  thunderstorm  centers  was  in  Africa.  The  round-the-world  path 
connecting  Africa  and  Massadhusetts  passes  close  to  the  extremities  of  the 
polar  caps  rather  than  squarely  through  them.  Therefore,  contrary  to  the 
21  February  PCA  discussed  above,  the  full  Impact  of  the  28  January  PCA 
was  probably  not  reflected  in  the  Massachusetts  measurement.  Roughly  half 
of  the  first  Fresnel  zone  of  the  propagation  path  followed  by  the  atmospheric 
noise  would  be  subjected  to  the  PCA-dlsturbed  ionosphere.  Field  (1978) 
showed  that,  in  such  a circumstance,  the  effect  of  the  PCA  would  be  about 
half  that  obtained  if  the  entire  Fresnel  zone  had  been  affected.  Thus, 
whereas  the  21  February  event  was  calculated  to  have  caused  a shift,  6f^, 
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of  between  -0.35  Hz  and  -0.6  Hz — depending  on  the  strength  of  the  event — 
the  corresponding  numbers  for  the  28  January  event  would  be  between  -0.18  Hz 
and  -0.3  Hz.  Again,  given  the  many  uncertainties,  these  results  agree 
tolerably  well  with  the  measured  value  (Table  2)  of  Af^  SS  -0.4  Hz. 

The  fact  that  the  data  indicates  a "moderate"  decrease  in  on  28 
January  versus  only  a "slight"  decrease  for  the  apparently  more  potent 
21  February  event  is  somewhat  puzzling,  although  not  too  significant  in 
view  of  the  qualitative  nature  of  the  data  on 

Effects  of  X-Ray  Flares 

Data  briefly  summarized  in  Sec.  Ill  show  minimal  correlation  between 
SIDs  and  the  frequencies  and  Q’s  of  Schumann  resonances.  Calculations 
using  the  solar-flare  model  ionospheres  discussed  above  give  very  slight 
differences  between  ambient  and  disturbed  values  of  S^^  and  0 at  Schumann- 
resonance  frequencies.  Thus,  in  a trivial  sense,  the  calculations  and  data 
are  consistent.  No  significant  effects  were  calculated  and  none  were 
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V.  CONCLUSIONS 


An  extensive  literature  survey  has  Identified  a number  of  low-frequency 
noise  measurements  relevant  to  the  effects  of  ionospheric  disturbances  on 
ELF  propagation.  These  measurements  pertain  to  the  behavior  of  earth- 
ionosphere  cavity  resonances  and  ELF/VLF  noise  intensity  at  the  time  of  1) 
the  Starfish  high-altitude  nuclear  detonation,  2)  more  than  thirty  solar 
flares,  and  3)  two  polar-cap-absorption  events  (PCAs).  Until  controlled- 
source  ELF  experiments  are  carried  out  under  disturbed  conditions,  these 
collateral  data  provide  the  only  means  by  which  codes  that  predict  ELF 
propagation  anomalies  in  nuclear  environments  can  be  compared  with 
experiment. 

A full-wave  propagation  code  similar  to  predictive  codes  developed 
under  DNA  auspices  has  been  used  above  to  compute  ELF  propagation  anomalies 
due  to  Starfish,  PCA  events,  and  a large  number  of  model  x-ray  flares.  The 
results  of  these  calculations  are  compared — to  the  extent  possible — to  the 
available  data.  Despite  the  inadequacies  of  the  data,  we  are  able  to 
verify  key  theoretical  predictions. 

Good  agreement  is  demonstrated  between  code  predictions  and  ELF/VLF 
noise  intensity  during  solar  flares.  Both  experiment  and  theory  consis- 
tently show  that  the  vast  majority  of  flares  cause  the  intensity  to  increase 
at  frequencies  between  30  Hz  and  1 kHz,  decrease  between  5 kHz  and  10  kHz, 
and  increase  between  10  kHz  and  30  kHz.  This  agreement,  which  occurs  over 
three  frequency  decades,  strongly  implies  the  correctness  of  the  predictive 
codes . 

Although  the  Schumann-resonance  band  (8-to-30  Hz)  is  somewhat  lower 
than  the  ELF  communications  band  (45-to-80  Hz),  both  bands  propagate  in  the 
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TEM  waveguide  mode.  Data  on  Schumann  resonances  may  therefore  be  meaning- 
fully compared  with  the  outputs  of  ELF  propagation  codes.  The  results  of 
such  a comparison  made  In  this  report  semlquantltatlvely  support  the 
validity  of  the  codes.  Specifically,  the  numerical  calculations  correctly 
predict  that  either  Starfish  or  the  PCAs  would  cause  the  peak  frequencies 
of  the  Schumann  resonances  to  shift  downward  by  0.5  Hz  to  1 Hz  and  the 
Q's  of  the  resonances  to  decrease.  The  shift  in  resonant  frequencies  Is 
due  to  the  depressed  Ionosphere  causing  the  phase  velocity  of  the  signal 
to  decrease,  and  the  decrease  In  Q Is  caused  by  an  Increase  In  the  attenua- 
tion rate  In  the  earth- Ionosphere  waveguide. 
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DYC,  J.  Kanrn 

Aerospace  Corp. 

Deputy  Chief  of  Staff 
Research,  Development,  & Acq. 
Department  of  the  Air  Force 
ATTN:  AFRDQ 


AHN:  F,  Morse 


Analytical  Systems  Engineering  Corp. 
ATTN:  Radio  Sciences 


Deputy  Chief  of  Staff 
Program  and  Analyses 
Department  of  the  Air  Force 
ATTN:  PACSC,  R.  Paul 


Berkeley  Research  Associates,  Inc. 
ATTN:  J.  Workman 


Electronic  Systems  Division,  AFSC 
ATTN:  XRW,  J.  Deas 
ATTN:  YSEA 


Boeing  Co. 
ATTN 
ATTN 
AHN 
AHN 
ATTN 


J.  Kenney 
G.  Keister 
D.  Murray 
G.  Hall 
S.  Tashird 


Foreign  Technology  Division,  AFSC 
ATTN:  TQTD,  B.  Ballard 
AHN:  Nils,  Library 


University  of  California  at  San  Diego 
ATTN:  H.  Booker 


Rome  Air  Development  Center,  AFSC 
AHN:  EEP 


Computer  Sciences  Corp. 
AHN:  H.  Blank 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

General  Electric  Co. 

Space  Division 

ATTN;  M.  Bortner 

General  Electric  Co. 

AHN;  F.  Reibert 

General  Electric  Company-TEMPO 
Center  for  Advanced  Studies 


ATTN 

ATTN 

ATTN 


W.  Knapp 
DAS  I AC 
B.  Gambill 


General  Research  Corp. 

Santa  Barbara  Division 
ATTN:  J.  Ise,  Jr. 

ATTN:  J.  Garbarino 

Geophysical  Institute 
University  of  Alaska 

ATTN;  Technical  Library 

GTE  Sylvania,  Inc. 

Electronics  Systems  Grp. -Eastern  Div. 
AHN;  M.  Cross 

Institute  for  Defense  Analyses 
ATTN:  J.  Aein 

International  Tel.  & Telegraph  Corp. 
ATTN:  Technical  Library 

JAYCOR 

ATTN:  S.  Goldman 

Johns  Hopkins  University 
Applied  Physics  Lab. 

ATTN;  T.  Potemra 
ATTN:  Document  Librarian 
ATTN;  T.  Evans 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Science  Applications,  Inc. 

ATTN:  D.  Sachs 


SRI  International 


ATTN 

ATTN 

ATTN 


G.  Carpenter 
0.  Neil  son 
G.  Price 


Lawrence  Livermore  Laboratory 
University  of  California 

ATTN:  Doc.  Con.  for  Technical  Information 
Dept.  Library 

M.l.T.  Lincoln  Lab. 

ATTN;  D.  Towle 

Mission  Research  Corp. 

ATTN:  D.  Sowle 
ATTN:  M.  Scheibe 

Mitre  Corp. 

ATTN:  C.  Callahan 
ATTN;  G.  Harding 

Mitre  Corp. 

ATTN:  M.  Horrocks 

Pacific-Sierra  Research  Corp. 

ATTN;  E.  Field,  Jr. 

R & D Associates 

ATTN:  R.  Lelevier 
ATTN:  B.  Gabbard 

Rand  Corp. 

ATTN:  C.  Crain 


